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ABSTRACT 

The preliminary design of multistage axial compressors in 
gas turbine engines is typically accomplished with mean-line 
methods. These methods, which rely on empirical correlations, 
estimate compressor performance well near the design point, 
but may become less reliable off-design. For land-based 
applications of gas turbine engines, off-design performance 
estimates are becoming increasingly important, as turbine plant 
operators desire peaking or load-following capabilities and 
hot-day operability. The current work develops a one- 
dimensional stage stacking procedure. This includes a newly- 
defined blockage term, which is used to estimate the off-design 
performance and operability range of a 13 -stage axial 
compressor. The new blockage term is defined to give 
mathematical closure on static pressure, and values of blockage 
are shown to collapse to a curve as functions of stage inlet flow 
coefficient and corrected speed. Utility of the stage stacking 
procedure is demonstrated by estimation of the minimum 
corrected speed which allows stable operation of the 
compressor. Further utility of the stage stacking procedure is 
demonstrated with a bleed sensitivity study, which estimates a 
bleed schedule to expand the compressor’s operating range. 

NOMENCLATURE 
Roman Symbols 

A an Geometrical annulus area, m 1 2 

A ef Effective flow area, m 2 

c p Specific heat at constant pressure, J/(kg K) 

k b Effective blockage factor 

M Mach number 

m Mass flow rate, kg/s 

N s Shaft rotational speed 


N c Corrected rotational speed 

P Total (stagnation) pressure, Pa 

P s] NASA standard sea-level pressure = 101,315 Pa 

p s Static pressure, Pa 

R Specific gas constant = 287.04 J/(kg K) 

r m Midspan radius, m 

T Total (stagnation) temperature, K 

T s i NASA standard sea-level temperature = 288.17 K 

U m Midspan rotor speed, m/s 

V z Axial flow velocity, m/s 

Greek Symbols 

a Absolute flow angle, deg 

y Ratio of specific heats 

// Adiabatic efficiency 

cp Flow coefficient 

iff Ideal work coefficient 

co Shaft angular frequency, rad/s 

Subscripts 

1 Stage inlet station 

2 Stage exit station 

INTRODUCTION 

Compressors in modern industrial gas turbine engines used 
for power generation are designed to operate continuously at a 
physical shaft speed corresponding to the utility frequency of 
either 50 Hz or 60 Hz. Industrial engines generally operate 
continuously for months while generating a fixed power output, 
and only require inspection and maintenance every few years. 
Typical maintenance intervals for these “baseload” engines are 
on the order of 25,000 hours or 800 starts [1]. Due to the 
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operation of these engines at a fixed shaft speed, stability issues 
arise. The engine power output requirement and shaft speed are 
fixed, whereas the compressor inlet temperature varies 
according to the local climate, causing a change in the corrected 
speed defined as 



( 1 ) 


Here N s is the physical shaft speed, T is the ambient, or 
compressor inlet total temperature, and T si is the NASA 
standard sea-level temperature of 288.17 K. It has been 
reported that industrial gas turbine engine power output can 
drop by as much as 0.9% for every 1 K rise in ambient 
temperature [2]. Care must be taken during the preliminary 
design of these compressors to ensure stable operation across 
the expected range of ambient temperatures. 

Besides the design requirements imposed by hot day 
operation, additional requirements are emerging due to the 
expected increase of renewable energy sources in the coming 
years [1]. Wind turbines and solar cells provide a varying 
amount of power dictated by local weather conditions. During 
cloudy days or days with low wind speed, gas turbine engines 
may be used to supplement power, especially during the hours 
when consumption demand, and thus the price of electricity, is 
high. This large demand generally exceeds the capacity that 
can be met by baseload engines, wind turbines, and solar cells. 
Future gas turbine engines must perform under a cyclic service 
duty to achieve this “peaking” capability. The future engine 
must be able to provide a range of power output at fixed shaft 
speed (as opposed to a fixed power output at fixed shaft speed) 
so that the engine output can follow demand throughout the 
day. This effectively requires the engine to be designed for 
multiple points of operation and presents an engineering 
challenge in the aerodynamic design of the compressor. 

In a preliminary aerodynamic design system, empirically 
based methods, referred to as mean-line or throughflow models, 
are used by engine manufacturers to establish an initial 
geometry for a new compressor [3-5]. These techniques make 
use of the breadth of experience and knowledge accrued from 
past compressor designs in the form of empirical or semi- 
empirical correlations to estimate parameters of the flow, such 
as blockage, loss, and deviation. These parameters have a large 
impact upon the performance and operability of the compressor. 
These mean-line or throughflow models can estimate the 
performance of a new compressor design quite well near its 
design point. However, these methods can become less reliable 
off-design, or for new designs which differ substantially from 
previous ones. Off-design performance estimates are becoming 
more important, as hot-day operation and load- folio wing or 
peaking capabilities become increasingly desirable for turbine 
plant operators. This poses a challenge during a preliminary 
design when critical performance estimates far from the design 
point are required. Thus, there exists a need for a methodology 


to generate credible estimates of off-design performance and 
operability of these compressors during preliminary design. 

With the increasing power and availability of computer 
resources, engine manufacturers can now use three-dimensional 
computational fluid dynamics (CFD) solvers to examine off- 
design operation of a compressor. These solvers rely less on 
empirical correlations than most mean-line and throughflow 
methods. A properly validated and calibrated CFD solver can 
thus provide insight into off-design compressor performance. 

The objective of this work is to develop a methodology 
which leads to credible estimates of aero -operability and aero- 
performance of axial multistage compressors at design and at 
off-design operating points for use in the preliminary design of 
a compressor whose design parameters are near that of an 
existing configuration. This objective is accomplished by 
developing a one-dimensional stage stacking procedure and 
using it to estimate design and off-design performance and 
operability limits of a 13 -stage axial flow compressor. This 
compressor is part of an industrial gas turbine engine used for 
power generation. The stage stacking procedure utilizes stage 
performance curves which are constructed using flow 
information extracted from a database of CFD simulations of 
the compressor. This database is generated using the APNAS A 
CFD code for multistage turbomachinery. 

DETAILS OF THE COMPRESSOR AND CFD MODEL 

The focus of the current work is a 13-stage axial flow 
compressor which is part of an industrial gas turbine engine. 
The design shaft speed for the compressor is 3600 RPM, 
corresponding to a 60 Hz utility frequency. The compressor 
consists of an inlet guide vane (IGV), 13 rotor-stator stages, and 
an outlet guide vane (OGV). All 28 blade rows are modeled in 
the current work. The IGV as well as stators 1, 2, and 3 
(referred to as variable guide vanes, or VGVs) have variable 
stagger angles. The remaining stators and OGV are 
cantilevered. 

A computational grid modeling the compressor was created 
using the mesh generator detailed by Mulac [6] . This generator 
produces an H-mesh for each blade row, with common axial 
and radial coordinates shared between the meshes. The meshes 
contain 51 nodes in the tangential direction (from the suction 
side of one blade to the pressure side of the next blade) and 5 1 
nodes in the radial direction (from hub to casing). Each blade 
is modeled with 61 axial nodes along the chord, from leading 
edge to trailing edge. The number of axial nodes from the inlet 
plane to the exit plane of the grid totals 3067 nodes. A 
meridional view of the compressor flow path is shown in 
Figure 1. 



FIGURE 1 : MERIDIONAL VIEW OF THE COMPRESSOR 
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The three VGVs were set to nominal settings, and the IGV 
was closed 2° from nominal in order to match the flow rate 
obtained from previous experimental tests of this compressor. 
These setting angles were fixed irrespective of simulated 
operating conditions. Stator platform leakages associated with 
the IGV and VGVs are not modeled in this work. This is 
because leakage estimates for this machine were extremely 
small compared to what is typical in compressors for aircraft 
engines. There are three casing bleeds in the compressor, 
located downstream of stators 5, 8, and 11. The removal of 
mass flow due to these bleeds is accounted for in the CFD code 
by using an aspiration model. Bleed flow was removed along 
an axisymmetric ring at the casing, downstream of the relevant 
stators. At design speed, the bleed flow rates are set to their 
nominal values. For part-speed operation, bleed rates are 
scaled by percent corrected speed. The physical geometry of 
the casing bleeds are not modeled in this work. 

CFD simulations of the compressor were generated using 
the APNASA turbomachinery code. The APNASA code, 
developed by Adamczyk, et al. [7], solves a set of equations 
known as the Average-Passage model, derived by Adamczyk 
[8], based on a sequence of mathematical averaging of the 
Navier-Stokes equations, continuity equation, energy equation, 
and the equation of state. The code employs a standard k-s 
turbulence model to close on the Reynolds stress term. Details 
of the code algorithm and governing equations can be found in 
[7-10]. 

An inlet boundary condition was prescribed at the first 
axial plane of the computational grid, at a location upstream of 
the IGV leading edge. This inlet boundary condition was 
specified as a radial profile of total pressure, total temperature, 
tangential flow angle, and radial flow angle. The inlet profile 
was obtained from an existing throughflow calculation for this 
machine, and was fixed irrespective of operating conditions. 
An exit boundary condition was prescribed at the last axial 
plane of the computational grid. This exit boundary condition 
was a static pressure specified at the hub using a simple radial 
equilibrium model to establish the spanwise flow distribution 
[7]. The simulation was throttled along a speed line by varying 
this exit hub static pressure. 

The corrected speed was varied by specifying the physical 
shaft speed, as opposed to changing the inlet temperature. This 
was done to emulate a series of surge tests that had previously 
been performed for this compressor. In these experiments, the 
IGV and VGV settings were fixed and mechanical shaft speed 
was reduced until the occurrence of a surge event. The 
reduction in speed was done at a rate such that clearances were 
minimally affected. It should be noted that this operation 
differs from routine operation of the engine, in which corrected 
speed varies according to changes in ambient temperature 
conditions while the physical shaft speed of the compressor is 
fixed according to the utility frequency of 60 Hz. The 
simulation results were compared to a limited set of stator 
leading edge total pressure and total temperature 
measurements, and agreement was judged to be satisfactory for 
the purposes of this work. 


STAGE STACKING PROCEDURE 

A one-dimensional procedure for estimating overall 
compressor performance, known as stage stacking, was 
described by Robbins and Dugan [11]. This method is 
frequently used in the preliminary design of axial compressors. 
The stage stacking procedure requires performance curves for 
every stage to be available. The performance curves consist of 
an ideal work coefficient and an adiabatic efficiency which are 
plotted as a function of stage inlet flow coefficient. The current 
work expands upon the procedure described by Robbins and 
Dugan [11] by defining and implementing a blockage term into 
the stage stacking procedure. If modifying an existing design, 
the performance curves may be generated from experimental 
tests. In the current work, the stage performance curves are 
generated using flow information extracted from CFD 
simulations, as will be detailed later. 

Stage Performance Parameters 

The stage ideal work coefficient if/ is defined as the ideal 
total enthalpy rise across the stage, non-dimensionalized by 
midspan wheel speed squared, as in 



( 2 ) 


In the current work, total temperature T / and total pressures 
Pi and P 2 are mass averaged values obtained from a CFD 
derived database. The subscripts 1 and 2 indicate the stage inlet 
and exit stations, respectively. The NASA standard sea-level 
total pressure and total temperature, P s! and T sh are 101,315 Pa 
and 288.17 K, respectively. The ratio of specific heats used 
here, y avg , is the arithmetic average of the ratios of specific heat 
at the stage inlet and stage exit stations, yj and y 2 , respectively. 
The specific heat at constant pressure c p is given by 


c 


p 


g_Yavg_ 
Yavg ~ 1 


( 3 ) 


Here, R is the specific gas constant for air, 287.04 J/(kg K). 
The midspan rotor speed at the stage inlet station U m i is given 
by 


U m , 1 = (4) 

Here, co is the angular frequency of the rotating shaft with 
units of rad/s and r m>1 is the midspan radius at the stage inlet 
station. Stage adiabatic efficiency rj is defined as the ratio of 
the ideal total enthalpy rise to the actual total enthalpy rise, as 
given by 
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Yavg 1 



The stage inlet flow coefficient cp 7 is defined as the ratio of 
corrected axial velocity to corrected midspan rotor speed, as 
given by 



The corrected axial velocity in Eq. (6) is solved from the 
following equation from [11]: 



The design value of stage inlet midspan absolute flow 
angle a 7 is used here, and m is the mass flow rate at the stage 
inlet station. The value of stage inlet mass flow rate is equal to 
the compressor inlet mass flow rate, less any applicable flow 
removal due to bleeds. The current work expands upon the 
method described in [11] regarding the area term A efl on the 
left-hand-side of Eq. (7). This term is generally taken to be the 
geometric annulus area, but more realistic values of corrected 
axial velocity are achieved if one accounts for flow blockage 
[11]. In the current work, the term A ef is taken as an effective 
flow area, defined such that it provides mathematical closure on 
static pressure p s>1 as given by the equations 



71 e/,1 
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—M 1 cos a ± 
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( 8 ) 
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Yi~ 1 
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(9) 


Combining the above equations and solving for A ef>1 gives 


A ef, 1 




Yi+l 



(Vs, l) 2 y x 

(EnY^—) 1 1 

1 2 

y j cos ui 

|UJ ; 

\yi - 1 


( 10 ) 


The static pressure term is taken as the area averaged 
value, extracted from CFD simulations of the compressor, but 
may alternatively be supplied by experimental tests or by 
correlations. Using this effective area term, an effective 
blockage term k b is defined as 


A e f 

k b = l--^ (11) 

™an 

Here A an is the geometric annulus area. Thus, the 
following information must be known to estimate compressor 
performance using the stage stacking procedure developed in 
the current work: 

1. Fully populated performance curves (iff, //, k b vs. <p 7 ) for 
each stage at the corrected speeds of interest. 

2. Geometrical annulus area at each stage inlet and exit 
station (A anJ and A an>2 ). 

3. Midspan radius at each stage inlet and exit station (r m 7 and 

^m, 2\ 

4. Design value of midspan absolute flow angle at each stage 
inlet and exit station (a 7 and a 2 ). 

Once this information is available, the stage stacking 
procedure is initiated by specifying the stage 1 inlet flow 
coefficient (p h inlet total pressure P h inlet total temperature T h 
and compressor corrected speed N c . Selection of these 
parameters sets the values of iff, //, k b>1 , and k b>2 for stage 1 
according to the stage 1 performance curves. With these 
parameters known, Eqs. (2-11) can be rearranged to solve for 
the stage 1 exit flow coefficient (p 2 , as detailed in [11, 12]. 

For a stage without bleed, such as stage 1, the exit station 
and the inlet station of the downstream stage defined at an 
identical location (i.e. midgap), and thus <p 2 of this stage is 
equal to the inlet flow coefficient of the following stage. This 
value of stage 2 inlet flow coefficient, along with the specified 
N c , sets the values of y, q, k bJ , and k h2 for stage 2 according to 
its performance curves. This process is repeated for each 
subsequent stage through the entire machine, yielding an 
estimate of the overall performance of the compressor. Note 
that for stages neighboring interstage bleeds, the exit station of 
the upstream stage is at a different location than the inlet station 
of the downstream stage. Thus, the value of cp 2 is corrected to 
account for mass flow removal due to the bleed and for any 
change in annulus area between the exit and inlet stations. 

Constructing Stage Performance Curves 

The stage stacking procedure defined above requires, 
among other information, the fully defined performance curves 
of every stage in the compressor. During preliminary design, 
the performance curves are generally estimated from empirical 
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correlations derived from test data. These correlations are well 
grounded near design point. However, at off-design conditions, 
their ability to estimate correct levels of performance is very 
often found to be lacking. In the current work, the stage 
performance curves were constructed using flow information 
extracted from APNASA CFD simulations. This database of 
CFD simulations had to be sufficiently extensive to allow a 
stage stacking method to be used to estimate compressor 
performance far from the design point. Simulations of the 
entire compressor will not, in general, yield the information 
required to fully populate the performance curves for all of its 
single stages. This is because certain stages will stall or choke, 
depending on operating conditions, well before others due to 
the way the stages are aerodynamically matched. In order to 
fully define the near- stall and near-choke portions of the 
characteristics for each stage, it is necessary to simulate small 
groups of stages in isolation. 

Since stage characteristics were constructed from 
simulations of small groups of stages, it was necessary to show 
that stage characteristics constructed from these simulations 
matched those constructed from simulations of the entire 
compressor. To this end, a design speed simulation of stages 1 
through 6 was generated and throttled along a speed line. Stage 
performance curves constructed from this simulation were 
compared to those constructed from a simulation of the entire 
compressor at design speed. The stage 1 performance curves of 
these two simulations are shown in Figure 2. 



FIGURE 2: STAGE 1 PERFORMANCE CURVES AT DESIGN 
SPEED - FULL COMPRESSOR VERSUS STAGES 1-6 


Points shown as open diamonds track the performance of 
stage 1 in the full compressor simulation and points shown as 
solid squares track the performance of stage 1 in a simulation of 
stages 1 through 6 only. The left-most points from both sets of 
data correspond to the near-stall points of the two simulations. 
Two major points are to be made from the results shown. 
Firstly, it is clear that the stage 1 characteristic constructed from 
the simulation of stages 1 through 6 lays on top of the stage 1 
characteristic constructed from the full compressor simulation. 


Therefore, it is valid to use results from simulations of groups 
of stages in place of full compressor simulations for the 
purposes of generating stage performance curves. Secondly, it 
is evident that simulations of the full compressor cannot 
adequately simulate the near- stall behavior of stage 1. At 
design speed, the full compressor simulation goes into 
numerical stall at a much greater flow coefficient than the 
simulation of stages 1 through 6. This is because a stage 
downstream of stage 6 is thought to be the stall limiting stage at 
design speed. Exclusion of this stage in the simulation of 
stages 1 through 6 allows stage 1 to be throttled to much lower 
flow coefficients as compared to the full compressor 
simulation. This demonstrates the need to simulate smaller 
groups of stages, rather than the full compressor, in order to 
fully populate the stage characteristics. 

Another example showing the limited flow range for a 
stage simulated within the full compressor is seen in Figure 3. 
This figure plots the performance of stage 10. Open diamonds 
show the performance of stage 10 within full compressor 
simulations, which are unable to establish the choking flow 
coefficient of stage 10. The choking flow coefficient was 
instead determined from simulations of stages 8 through 10 in 
isolation, as shown by the solid squares in Figure 3. Choking 
of the stage is identified by the near-vertical slope of the 
efficiency curve at high flow coefficient. 



FIGURE 3: STAGE 10 PERFORMANCE CURVES AT DESIGN 
SPEED - FULL COMPRESSOR VERSUS STAGES 8-10 

In simulating a group of embedded stages (e.g. stage 10) in 
this manner, it became necessary to include at least two stages 
upstream of the stage of interest. In [12], it was shown that the 
two upstream stages act as flow generators, which establish the 
multistage environment for the stage of interest. This is 
consistent with industry practices. Two stages upstream were 
adequate for generating the correct inlet conditions to the stage 
of interest as the simulation was throttled. Additionally, some 
simulations included at least one stage downstream, but it was 
found that in many cases, the downstream stage could be 
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removed without impacting the characteristic curves, as shown 
in [12]. 

The stage performance curves generated in this manner 
displayed results consistent with [11]. As noted in [11], stages 
with inlet relative Mach number greater than about 0.75 have 
discrete curves, iff and tj as a function of <pi, at different 
corrected speeds. In the current compressor, stages 1 through 4 
have inlet Mach number greater than 0.75. Figure 4 shows the 
performance curves for stage 1, which exhibit this dependence 
on corrected speed and on c 


c 

| 1.20 
it 

Q) 

5 115 
| 1.10 

1 105 
I 1.00 

I 0 - 9 * 



70 0.72 0.74 0.76 0.78 0.80 0.82 0.84 


Normalized Stage Inlet Flow Coefficient 


FIGURE 4: STAGE 1 PERFORMANCE CURVES AT VARIOUS 
CORRECTED SPEEDS 


Furthermore, as noted in [11], stages with inlet relative 
Mach number less than 0.75 have performance curves which 
collapse to a single line irrespective of corrected speed. This 
was observed for stages 5 through 13, which have inlet relative 
Mach numbers less than 0.75. The collapse of these curves is 
shown in Figure 5, which shows the performance curves of 
stage 5 . 


It is clear from Figure 5 that the performance curves for 
this stage have negligible dependence on corrected speed. 
Thus, performance curves of stages 1 through 4 were generated 
at various corrected speeds, whereas the performance curves of 
stages 5 through 13 could be fully populated with only design 
speed simulations. 

Collapse of Blockage Curves 

A noteworthy result observed when constructing the stage 
performance curves was the behavior of the blockage parameter 
k b . When plotted as a function of (p x this parameter, like if/ and 
/ 7 , showed a dependence on corrected speed for those stages 
with inlet relative Mach number greater than 0.75 (i.e. stages 1 
through 4). This is shown in Figure 6, which plots the effective 
blockage at the exit of stage 1 k b>2 as a function of stage 1 inlet 
flow coefficient (pi for various corrected speeds. It should be 
noted that the blockages shown here have been normalized by 
the value of stage 1 blockage near the design point. 



FIGURE 6: STAGE 1 BLOCKAGE CURVES AT VARIOUS 
CORRECTED SPEEDS 



FIGURE 5: STAGE 5 PERFORMANCE CURVES AT VARIOUS 
CORRECTED SPEEDS 


Furthermore, these blockage curves collapsed to a single 
line for those stages with inlet relative Mach number less than 
0.75 (i.e. stages 5 through 13) irrespective of corrected speed. 
This is shown in Figure 7, which shows the stage 5 exit 
blockage curves at various corrected speeds. 
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FIGURE 7: STAGE 5 BLOCKAGE CURVES AT VARIOUS 
CORRECTED SPEEDS 


This correlation of blockage with flow coefficient and 
corrected speed is a result which, to the authors’ knowledge, 
has not been published in the open literature. This correlation 
allowed for the implementation of the blockage curves into the 
stage stacking procedure described in [11], whereby the 
selection of inlet flow coefficient provided values of effective 
blockage in addition to values of the ideal work coefficient and 
adiabatic efficiency of a stage. 

A heuristic argument is made here to understand the 
collapse of these blockage curves. As shown in [12], the 
blockage parameter was shown to be consistent with an existing 
blockage metric defined by Khalid, et al. [13] over a wide range 
of operating conditions. The blockage parameter defined in 
[13] is a measure of regions of flow where the gradient of a 
velocity component is large. This is the same region of flow 
where the dissipation function, and thus the entropy production, 
is large. It can be shown from the equation of efficiency, Eq. 
(5), that the entropy production for stages with Mach number 
less than 0.75 is only a function of flow coefficient. Thus, one 
may conclude that the collapse of the blockage curves as shown 
in the present work is a generic result. 

APPLICATIONS OF STAGE STACKING 

The typical use of stage stacking methods is to generate 
stage and overall performance estimates for a given compressor 
at specified operating conditions. In the current work, the stage 
stacking procedure was used to estimate off-design compressor 
performance at 97% speed, as detailed in [12]. These overall 
compressor performance estimates yielded total pressure ratio 
to within 1.6% and adiabatic efficiency within 0.3 points of 
CFD simulations of the entire compressor. The comparison of 
CFD and stage stacking results was done as a plausibility check 
on the stage stacking procedure. 

In addition, the stage stacking procedure was used to 
estimate the minimum corrected speed for stable compressor 


operation with fixed IGV/VGV angles and fixed bleed rates. 
The details of this study follow. 

Estimation of Minimum Stable Corrected Speed 

The stage stacking procedure was used to estimate the 
minimum corrected speed which allowed stable compressor 
operation with geometry and bleed rates fixed. This was done 
by performing the stage-by-stage calculation, described 
previously, and picking near- stall values of stage 1 inlet flow 
coefficient along the various performance curves at different 
corrected speeds. The goal of the stage stacking calculations 
was to find the corrected speed which would result in a choking 
flow coefficient for the last stage in the compressor. Choked 
flow in the rear stage is typical of a compressor operating in 
isolation near minimum speed. After iteratively reducing 
corrected speed in the stage stacking calculations, it was found 
that operation at part-speed “N%” was characterized by flow 
coefficients near stall at the compressor inlet, and near choke at 
the compressor exit. Part- speed N% was found to be a 
corrected speed significantly less than design speed. 



FIGURE 8: STAGE 1 PERFORMANCE CURVES WITH STAGE 
STACKING ESTIMATE AT N% SPEED 

Figure 8 shows the stage 1 N% performance points as open 
squares. The solid diamond depicts the stage stacking 
estimated performance for the stage. As shown, the mi nimum 
(near-stall) value of flow coefficient was specified for the stage 
1 inlet flow coefficient in this stage stacking calculation. 
Selection of this flow coefficient determined the ideal work 
coefficient, adiabatic efficiency, and effective blockage 
according to the stage performance 1 curves. The stage-by- 
stage calculation was performed, resulting in a near-choke 
value of stage 13 inlet flow coefficient, as shown in Figure 9. 
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FIGURE 9: STAGE 13 PERFORMANCE CURVES WITH 
STAGE STACKING ESTIMATE AT N% SPEED 


A separate stage stacking calculation was performed for 
part-speed “N2%”, which is a corrected speed 2% less than 
part- speed N%. A near-stall value of stage 1 inlet flow 
coefficient was picked along the N2% speed stage 1 
performance curves to initiate the stage stacking calculation, as 
shown in Figure 10. 



FIGURE 10: STAGE 1 PERFORMANCE CURVES WITH 
STAGE STACKING ESTIMATE AT N2% SPEED 


As the stage-by-stage calculation progressed through the 
machine, it was noted that the stage stacking estimate of stage 9 
performance indicated that this stage was nearly choked, as 
shown in Figure 1 1 . This is readily seen on the efficiency plot, 
which shows that the stage stacking point is on the near- vertical 
portion of the curve. 



FIGURE 11 : STAGE 9 PERFORMANCE CURVES WITH 
STAGE STACKING ESTIMATE AT N2% SPEED 


Attempting to continue the stage stacking procedure 
beyond stage 9 proved to be futile. This was because the stage 
10 inlet flow coefficient estimated by the stage stacking 
procedure was greater than the maximum possible value 
attainable by the stage, as shown in Figure 12. 



FIGURE 12: STAGE 10 PERFORMANCE CURVES WITH 
STAGE STACKING ESTIMATE AT N2% SPEED 


The stage 10 inlet flow coefficient estimated via stage 
stacking was far greater than the maximum, or choking, flow 
coefficient. Stage 10 cannot operate at such a high flow 
coefficient. In order to achieve operability at part-speed N2%, 
the stage 10 inlet flow coefficient must decrease. For fixed 
geometry and bleed rates, this is only possible if the compressor 
inlet flow coefficient is reduced. However, a near-stall value of 
stage 1 inlet flow coefficient was specified in this case. The 
compressor would likely stall if it was to operate at a lower 
inlet flow coefficient. Hence at part-speed N2%, there is no 
operating range for this compressor for bleed rates and 
IGV/VGV settings fixed to their nominal values. Thus, the 
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minimum stable speed of the compressor is estimated to be 
greater than N2%, and equal to or less than N% speed. The 
final example, which follows, details a bleed sensitivity study. 

Operability and Bleed Sensitivity Study 

A bleed sensitivity study was conducted to estimate the 
stage 5 bleed rate required such that stage 13 operated near- 
choke, and all other stages operated within their stable 
operating ranges at part- speed N2%. This was accomplished by 
iteratively varying stage 5 bleed rate within the stage stacking 
procedure and performing the stage-by-stage calculation with 
the stage 1 inlet flow coefficient fixed to its near- stall value for 
N2% speed. These iterations were performed until the stage 
stacking procedure estimated a near-choke flow coefficient for 
stage 13, with all upstream stages operating in their stable 
operating ranges. The resulting stage 5 bleed rate was 
approximately 7 times greater than the nominal rate. 

The value of exit static pressure resulting from the stage 
stacking estimate was specified as the exit static pressure for a 
CFD simulation at N2% speed using the new value of stage 5 
bleed rate, as detailed in [12]. Stage performance parameters 
calculated from this CFD simulation were then compared to the 
stage stacking estimate. Figure 13 shows the stage 1 
performance curves. 


Stage 1 Characteristic, N2% Speed 
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FIGURE 13: STAGE 1 PERFORMANCE CURVES WITH 
STAGE STACKING AND CFD RESULTS AT N2% SPEED 


The near-stall value of flow coefficient is denoted by the solid 
diamond. The CFD simulation is shown as a solid triangle in 
Figure 13, and it indicates that the CFD simulation converged 
to a slightly lower stage 1 inlet flow coefficient at N2% speed 
as compared to the stage stacking estimate. The result of 
increasing the stage 5 bleed rate for stage 9 is shown in Figure 
14. 



FIGURE 14: STAGE 9 PERFORMANCE CURVES WITH 
STAGE STACKING AND CFD RESULTS AT N2% SPEED 


The stage stacking estimate of stage 9 inlet flow coefficient 
with nominal stage 5 bleed rate (shown as a solid circle) is 
compared to that with increased stage 5 bleed rate (shown as a 
solid diamond). The increase in stage 5 bleed rate resulted in a 
significantly reduced estimate of stage 9 inlet flow coefficient. 
Examining the stage 13 performance curves in Figure 15 shows 
that the CFD simulation resulted in stage 13 being choked, 
which was the goal of the bleed sensitivity study. 
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FIGURE 15: STAGE 13 PERFORMANCE CURVES WITH 
STAGE STACKING AND CFD RESULTS AT N2% SPEED 


Thus, the stage stacking procedure is able to describe a 
bleed schedule which could extend the operating range of this 
compressor down to part-speed N2%, which is a corrected 
speed 2% lower than N% speed. Increasing the stage 5 bleed 
rate in the compressor acted to reduce the flow coefficient for 
choked stages within the middle of the machine. The stage 
stacking calculation estimated the compressor inlet corrected 
flow to within 0.3% of the resulting CFD simulation. Overall 
total pressure ratio was estimated to within 1.3% of the CFD 
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result, and overall adiabatic efficiency was estimated to within 
0.42 points of the CFD result. 

As a final step in this operability study, the stage 5 bleed 
rate in the CFD simulation was reduced from 7 times the 
nominal value back down to the nominal value. This resulted 
in numerical stall in the front stages of the compressor. This 
result corroborated the estimated minimum stable corrected 
speed from the stage stacking calculations, as there was no 
operating range at part- speed N2% with nominal IGV/VGV and 
bleed settings. 

CONCLUSIONS 

The major conclusions of this research effort can be 
summarized as follows: 

• A modified stage stacking procedure was developed which 
attempted to account for flow blockage through a novel 
mathematical closure on static pressure to calculate 
effective flow area. 

• The blockage term, like the ideal work coefficient and 
adiabatic efficiency, was found to be a function of: 

o inlet flow coefficient and corrected speed, for those 
stages with inlet relative Mach number greater than 0.75 
(stages 1-4). 

o inlet flow coefficient only, for those stages with inlet 
relative Mach number less than 0.75 (stages 5-13). 

• To fully define the stage performance curves at the near- 
stall and near-choke sides using CFD results, groups of 
stages must be simulated in isolation. The complete flow 
range of a given stage cannot, in general, be modeled in 
simulations of the full compressor. 

• The stage stacking procedure yielded: 

o good off-design estimates of overall compressor 
performance as compared to CFD simulations, 
o a good estimate for the minimum corrected speed for 
stable operation as compared to experiment, 
o a good estimate of a bleed schedule to expand the part- 
speed operability of the compressor as compared to CFD 
simulations. 

The primary use of the proposed stage stacking procedure 
would be during the design of a compressor which implements 
incremental changes on an existing design; for example, when 
the designer is interested in the effect on overall compressor 
performance after changing bleed flow rates, IGV setting 
angles, or blade shape of a single stage. 
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